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Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for revisions 
in the picture of cosmic ray transport. In this paper we use recently obtained scaling laws for MHD 
modes to obtain the scattering frequency for cosmic rays. We account for the turbulence cutoff 
arising from both collisional and collisionless damping. We obtain the scattering rate and show that 
fast modes provide the dominant contribution to cosmic ray scattering for the typical interstellar 
conditions in spite of the fact that fast modes are subjected to damping. We determine how the 
efficiency of the scattering depends on the characteristics of ionized media, e.g. plasma /3. We show 
that streaming instability is suppressed by the ambient MHD turbulence. 



INTRODUCTION 



The propagation of cosmic rays (CRs) is affected 
by their interaction with magnetic field. This field is 
turbulent and therefore, the resonant interaction of 
cosmic rays with MHD turbulence has been discussed 
by many authors as the principal mechanism to scat- 
ter and isotropize cosmic rays. Although cosmic ray 
diffusion can happen while cosmic rays follow wander- 
ing magnetic fields (|14|). the acceleration of cosmic 
rays requires efficient scattering. For instance, scat- 
tering of cosmic rays back into the shock is a vital 
component of the first order Fermi acceleration. 



While most investigations are restricted to Alfven 
modes propagating along an external magnetic field 
(the so-called slab model of Alfvenic turbulence), 
obliquely propagating MHD modes have been in- 
cluded in |17| and later studies (25|. A more com- 
plex models were obtained by combining the results of 
the Reduced MHD with parallel slab-like modes have 
been also considered |3j . Here we attempt to use mod- 
els that are motivated by the recent studies of MHD 
turbulence (0], see f° r a review and references 



therein) . 



The efficiency of scattering depends on turbulence 
anisotropy [2(j • Therefore the calculations of CR scat- 
tering must be done using a realistic MHD turbu- 
lence model. An important attempt in this direction 
was carried out in |fj. However, only incompressible 
motions were considered. On the contrary, ISM is 
highly compressible. Compressible MHD turbulence 
has been studied recently (see review by [10| and ref- 
erences therein) . j2£| addressed the scattering by fast 
modes. But they did not consider the damping, which 
is essential for fast modes. In this paper we include 
various damping processes which can affect the fast 
modes. 



II. MHD TURBULENCE AND ITS DAMPING 
A. Model of MHD turbulence 

Analogous to MHD perturbations that can be de- 
composed into Alfvenic, slow and fast waves with well- 
defined dispersion relations, MHD perturbations that 
characterize turbulence can be separated into distinct 
modes. The separation into Alfven and pseudo- Alfven 
modes, which are the incompressible limit of slow 
modes, is an essential element of the GS95 0] model. 
Even in a compressible medium, MHD turbulence is 
not an inseparable mess in spite of the fact that MHD 
turbulence is a highly non-linear phenomenon [9j, |21| . 
The actual decomposition of MHD turbulence into 
Alfven, slow and fast modes was addressed in [jj , 
who used a statistical procedure of decomposition in 
the Fourier space, where the basis of the Alfven, slow 
and fast perturbations was defined. 

Unlike hydrodynamic turbulence, Alfvenic one is 
anisotropic, with eddies elongated along the magnetic 
field. On the intuitive level it can be explained as 
the result of the following fact: it is easier to mix 
the magnetic field lines perpendicular to the direction 
of the magnetic field rather than to bend them. How- 
ever, one cannot do mixing in the perpendicular direc- 
tion to very small scales without affecting the parallel 
scales. This is probably the major difference between 
the adopted model of Alfvenic perturbations and the 
Reduced MHD A\ . In the GS95 model as well as in its 
generalizations for compressible medium mixing mo- 
tions induce the reductions of the scales of the parallel 
perturbations. 

The corresponding scaling can be easily obtained. 
For instance, calculations in prove that motions 
perpendicular to magnetic field lines are essentially 
hydrodynamic. As the result, energy transfer rate due 
to those motions is constant Ek ~ v\jT^, where is 
the energy eddy turnover time ~ (vkk±) , where k± 
is the perpendicular component of the wave vector 
k. The mixing motions couple to the wave-like mo- 
tions parallel to magnetic field giving a critical balance 
condition, i.e., k±Vk ~ feii V^, where kit is the paral- 
lel component of the wave vector k, Va is the Alfven 
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speed 0. From these arguments, the scale dependent 

2 /3 

anisotropy fcii oc fc ± and a Kolmogorov-like spectrum 
for the perpendicular motions Vk oc fc -1 / 3 can be ob- 
tained jlflj . 

It was conjectured in [2l| that GS95 scaling should 
be approximately true for Alfven and slow modes in 
moderately compressible plasma. For magnetically 
dominated, the so-called low (3 plasma, |!| showed 
that the coupling of Alfvenic and compressible modes 
is weak and that the Alfvenic and slow modes follow 
the GS95 spectrum. This is consistent with the anal- 
ysis of HI velocity statistics |2j| as well as with 
the electron density statistics U Calculations in 01 
demonstrated that fast modes are marginally affected 
by Alfven modes and follow acoustic cascade in both 
high and low (3 medium. In what follows, we consider 
both Alfven modes and compressible modes and use 
the description of those modes obtained in 0, ^} to 
study CR scattering by MHD turbulence. 



B. Damping of Turbulence 

In many earlier papers Alfvenic turbulence was con- 
sidered by many authors as the default model of inter- 
stellar magnetic turbulence. This was partially moti- 
vated by the fact that unlike compressible modes, the 
Alfven ones are essentially free of damping in fully 
ionized medium 35]. However, it was shown that com- 
pressible fast modes are particularly important for 
cosmic ray scattering [3ll l33j . For them damping is 
essential. 

At small scales turbulence spectrum is altered by 
damping. Various processes can damp the MHD mo- 
tions (sec 33] for details). In partially ionized plasma, 
the ion-neutral collisions are the dominant damp- 
ing process. In fully ionized plasma, there are basi- 
cally two kinds of damping: collisional or collisionless 
damping. Their relative importance depends on the 
mean free path in the medium. If the wavelength is 
larger than the mean free path, viscous damping dom- 
inates. If, on the other hand, the wavelength is smaller 
than mean free path, then the plasma is in the colli- 
sionless regime and collisionless damping is dominant. 

To obtain the truncation scale, the damping time 
r^ 1 should be compared to the cascading time Tfc. 
As we mentioned earlier, the Alfvenic turbulence cas- 
cades over one eddy turn over time (k^Vk)^ 1 ~ 
(fen Va) — 1 . The cascade of fast modes takes a bit 
longer: 



r fc = u/k 2 v 2 = (k/L)- 1 / 2 x V ph /V 2 , 

where V is the turbulence velocity at the injection 
scale, V p h is is the phase speed of fast modes and 
equal to Alfven and sound velocity for high and low (3 
plasma, respectively . If the damping is faster than 
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FIG. 1: Damping scale vs. the angle 9 between k and B 
in halo, HIM and WIM. The peak on the curve of HIM 
(dashdot line) is smeared out by randomization of both k 
and B 13^1 . 



the turbulence cascade, the turbulence gets truncated. 
Otherwise, for the sake of simplicity, we ignore the 
damping and assume that the turbulence cascade is 
unaffected. As the transfer of energy between Alfven, 
slow and fast modes of MHD turbulence is suppressed 
at the scales less than the injection scale, we consider 
different components of MHD cascade independently. 

We get the cutoff scale k c by equating the damping 
rate and cascading rate r/Td — F Then we check 
whether it is self-consistent by comparing the k c with 
the relevant scales, e.g., injection scale, mean free path 
and the ion gyro-scale. 

Damping is, in general, anisotropic, i.e., the linear 
damping (see [3^] Appendix A) depends on the angle 
9 between the wave vector k and local direction of 
magnetic field B. Unless randomization of 9 is com- 
parable to the cascading rate the damping scale gets 
angle-dependent [3^. The angle 9 varies because of 
both the randomization of wave vector k and the wan- 
dering of magnetic field lines ( see [33 for detail). With 
this input at hand, one can determine the turbulence 
damping scales given a medium. 



III. INTERACTIONS BETWEEN 
TURBULENCE AND PARTICLES 

Basically there are two types of resonant interac- 
tions: gyroresonance acceleration and transit accel- 
eration (henceforth TTD). The resonant condition is 
u) — k\\v/Jt — nQ (n = 0,±1,2...), where u> is the wave 
frequency, CI — Qq/j is the gyrofrequency of relativis- 
tic particle, /1 = cos£, where £ is the pitch angle of 
particles. TTD formally corresponds to n = and it 
requires compressible perturbations. 

We employ quasi-linear theory (QLT) to obtain our 
estimates. If mean magnetic field is larger than the 
fluctuations at the injection scale, we may say that the 
QFT treatment we employ defines parallel diffusion. 
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Taking into account only the dominant interaction at 
n = ±1. we obtain the Fokker-Planck coefficients (see 

also Him), 



our results orders of magnitude larger than those 38] 
in [£j for the most of energies considered. However, 
the scattering frequency, 



m 2 Vl 



M nn {k) + M cc (k) 
K nn (k) + K cc (k) 



fi 



fi 

. k±v± 
fi 





' M TC£ (k) " 


+ e^ 2 * 


' Af £TC (k) " 




} 








. K cn (k) _ 







,(1) 



where £ = 1 for Alfven modes and £ = for fast 
modes, k min = fc mQ2 , = fi /«th corresponds to 

the dissipation scale, m = jrriH is the relativistic mass 
of the proton, v±_ is the particle's velocity component 
perpendicular to Bo, 4> — arctan(fc J/ / k x ) , C,TZ — (x ± 
iy) I y/2 represent left and right hand polarization. 

The delta function <5(fc||^|| — u + nfl) approximation 
to real interaction is true when magnetic perturba- 
tions can be considered static[37j ■ For cosmic rays, 
fc||i»H 3> <jJ so that the resonant condition is essentially 
ktiVfi — nfl = 0. From this resonance condition, we 
know that the most important interaction occurs at 



Q/vu. This is generally true except for 



small (j, (or scattering near 90°) 



2£W(1 - /i 2 



(3) 



are still much smaller than the estimates for isotropic 
and slab model j^Il- As the anisotropy of the Alfven 
modes is increasing with the decrease of scales, the 
interaction with Alfven modes becomes more efficient 
for higher energy cosmic rays. When the Larmor ra- 
dius of the particle becomes comparable to the injec- 
tion scale, which is likely to be true in the shock region 
as well as for very high energy cosmic rays in diffuse 
ISM, Alfven modes get important. 

It's worthwhile to mention the imbalanced cascade 
of Alfven modes Q . Our basic assumption above was 
that Alfven modes were balanced, meaning that the 
energy of modes propagating one way was equal to 
that in opposite direction. In reality, many turbulence 
sources are localized so that the modes leaving the 
sources are more energetic than those coming toward 
the sources. The energy transfer in the imbalanced 
cascade occurs at a slower rate, and the Alfven modes 
behave more like waves. The scattering by these im- 
balanced Alfven modes could be more efficient. How- 
ever, as the actual degree of anisotropy of imbalanced 
cascade is currently uncertain, and the process will be 
discussed elsewhere 391. 



IV. SCATTERING OF COSMIC RAYS 

A. Scattering by Alfvenic turbulence 

As we discussed above, Alfven modes are 
anisotropic, eddies are elongated along the magnetic 
field, i.e., k± > k\\. The scattering of CRs by Alfven 
modes is suppressed first because most turbulent en- 
ergy goes to fcj_ due to the anisotropy of the Alfvenic 
turbulence so that there is much less energy left in 



the resonance point ku 



Further- 



more, k± 3> fey means k± 3> r^ 1 so that cosmic ray 
particles have to be interacting with lots of eddies in 
one gyro period. This random walk substantially de- 
creases the scattering efficiency. The scattering by re- 
alistic Alfvenic turbulence was studied in [3J • In case 
that the pitch angle £ not close to 0, the analytical 
result is 



'pp 



Q-gT[6.5, k m axk\\ , res L ' 



m 2 Vl 



(2) 



where T[a, z] is the incomplete gamma function. The 
presence of this gamma function in our solution makes 



B. Scattering by fast modes 

The contribution from slow modes is no more than 
that by Alfven modes since the slow modes have 
the similar anisotropies and scalings. More promis- 
ing are fast modes, which are isotropic. With fast 
modes there can be both gyroresonance and transit- 
time damping (TTD). 

TTD happens due to the resonant interaction with 
parallel magnetic mirror force. The advantage of TTD 
is that it doesn't have a distinct resonant scale associ- 
ated with it. TTD is thus an alternative to scattering 
of low energy CRs which Larmor radii are below the 
damping scale of the fast modes. Moreover, we shall 
show later that TTD can contribute substantially to 
cosmic ray acceleration (also known as the second or- 
der Fermi acceleration). This can be crucial in some 
circumstances, e.g., for 7 ray burst (0), and acceler- 
ation of charged particles [32] • Different from gyrores- 
onance, the resonance function of TTD is broadened 
even for CRs with small pitch angles. The formal reso- 
nance peak fcy/fe = Vph/vu favors quasi-perpendicular 
modes. However, these quasi-perpendicular modes 
cannot form an effective mirror to confine CRs be- 
cause the gradient of magnetic perturbations along 
the mean field direction V11B is small. Thus the reso- 



2527 



4 



22nd Texas Symposium on relativistic Astrophysics at Stanford University, Dec. 13-17, 2004 



gyroresonance 
gyro analytical 



gyroresonance 

gyro analytical 

- - TTD 











10° E k (GeV) 10 2 " 10 - 2 10° E k (GeV) ,„* 



TTD 
TTD 

gyroresonance 
gyro analytical 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



FIG. 2: Scattering frequency v given by Eq.J^J vs. the ki- 
netic energy Ek of cosmic rays (a) in halo, (b) in WIM. The 
'x' lines refer to scattering by gyroresonance and solid lines 
are the corresponding analytical results given by Eq.Q. 
The dashed line are the results for TTD (from 



FIG. 3: The scattering frequency v vs. cosine of pitch 
angle fi of lGeV CR (a) in halo, (b) in WIM. The 'x' lines 
refer to scattering by gyroresonance and solid lines are 
the corresponding analytical results given by Eq.Q. The 
dashed and diamond lines are the results for TTD adopting 
5 function and Breit-Wigner function respectively (from 



nance peak is weighted out and the Breit-Wigner-type 
[271 resonance function should be adopted. 

Here we apply our analysis to the various phases 
of ISM. In low (3 medium, both collisionless and col- 
lisional damping increases with 8 unless 6 is close to 
7r/2. This means that those quasi-parallel fast modes 
are least damped. For these modes the argument for 
the Bessel function in Eq.JTJ is k± tan£/fc|i res < 1 
unless £ is close to 90°. So we can take advantage 
of the anisotropy of the damped fast modes and use 
the asymptotic form of Bessel function for small argu- 
ment J n {x) ~ (x/2) n /n\ to obtain the corresponding 
analytical result for this case: 



D 



7r(»^)°- 5 (l-/i 2 ) 
l)"i)/7 



(4) 



In high f3 medium and partially ionized medium, 
fast modes are subjected to severe damping and trun- 
cated at scale larger than the resonant scale of CRs, 
~ Q/c. The only available scattering mechanism is 
TTD, which can still provide much more efficient scat- 
tering for CRs comparing with Alfven modes (see (3j| 
for detail). 

A special case is that the cosmic rays propagate 
nearly perpendicular to the magnetic field, so called 
the 90° scattering problem. It should be noted that 
with resonance broadening associated with the tur- 
bulence, the requirement vn — V p h/cos9 > V p h is 
relieved. The contribution from TTD thus becomes 
dominant for sufficiently small fi. However, since 
quasi-linear approximation is not accurate in this 
regime, proper calculation should be carried out with 
non-linear effects taken into account 0, |2^| . 



V. COSMIC RAY SELF CONFINEMENT BY 
STREAMING INSTABILITY 

When cosmic rays stream at a velocity much larger 
than Alfven velocity, they can excite by gyroresonance 
MHD modes which in turn scatter cosmic rays back, 
thus increasing the amplitude of the resonant mode. 
This runaway process is known as streaming insta- 
bility. It was claimed that the instability could pro- 
vide confinement for cosmic rays with energy less than 
~ 10 2 GeV p|- However, this was calculated in an 
ideal regime, namely, there was no background MHD 
turbulence. In other words, it was thought that the 
self-excited modes would not be appreciably damped 
in fully ionized gas. 

This is not true for turbulent medium, however. 
pjlj pointed out that the streaming instability is 
partially suppressed in the presence of background 
turbulence[4(|. More recently, detailed calculations 
of the streaming instability in the presence of back- 
ground Alfvenic turbulence were presented in \H\ . 

For interaction with fast modes, it happens at the 
rateT fc ~ {k/ L)- l / 2 \UJV 2 (see Ea. lHTBl 'l. By equat- 
ing the growth rate |23 , 



tV/x o ^i— t 1 , ^stream - 

r(fc) = n (-1 + — — , 

n„ Va 



(5) 



and the damping rate Ea. (|IIB|> . we can find that the 
streaming instability is only applicable for particles 
with energy less than 

lmax ~ 1.5 x 10- 9 [n-\V ph /V)(Lvn /V 2 ) - 5 }^-\ 

(6) 

which for HIM, provides ~ 20GeV if taking the injec- 
tion speed to be V ~ 25km/s. 

One of the most vital cases for streaming instabil- 
ity is that of cosmic ray acceleration in strong MHD 
shocks. Such shocks produced by supernovae explo- 
sions scatter cosmic rays by the postshock turbulence 
and by preshock magnetic perturbations created by 
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cosmic ray streaming. The perturbations of the mag- 
netic field may be substantially larger than the regular 
magnetic field. The corresponding nonlinear growth 
rate is ([H): 

T(k) 

1 ' cV A ( 1?nax - (1 + a)-i)(fcr 90 )°(l + AU {1 ~ a)/2 ' 

(7) 

where e is the ratio of the pressure of CRs at the shock 
and the upstream momentum flux entering the shock 
front, U is the shock front speed, a — 4 is the spectrum 
index of CRs at the shock front, r g o — c/Q,q, A to t is the 
dimensionless amplitude of the random field Af ot = 
SB 2 /Bl 

Magnetic field itself is likely to be amplified through 
an inverse cascade of magnetic energy at which per- 
turbations created at a particular k diffuse in k space 
to smaller k thus inducing inverse cascade. As the re- 
sult the magnetic perturbations at smaller k get larger 
than the regular field. As the result, even if the insta- 
bility is suppressed for the growth rate given by eq. 
(0 it gets efficient due to the increase of perturbations 
of magnetic field stemming from the inverse cascade. 

Whether or not the streaming instability is efficient 
in scattering accelerated cosmic rays back depends on 
whether the growth rate of the streaming instability is 
larger or smaller than the damping rate. The precise 
picture of the process depends on yet not completely 
clear details of the inverse cascade of magnetic field. 
If, however, we assume that the small scale driving 
provides at the scales of interest isotropic turbulence 
the nonlinear damping happens on the scale of one 
eddy turnover time. Assuming the shock front speed 
U is low, we attain the maximum energy of particles 
accelerated in the shock by equating the growth rate 
of the instability Eq.{7|) to the damping rate due to 
turbulence cascade Ea. ljll B|) : 

^ MLeBpf^P /(0 . 5+a) 

_ ^ m 0.5y2 c 2 > ■ 

From this we can estimate jmax — 2 x 10 7 (t/kyr)~ 9 / 4 
in HIM. 

As shown in the Eq.© and JJJ, the growth rate 
depends on the CRs' density. In those regions where 



high energy particles are produced, e.g., shock fronts 
in ISM, 7 ray burst, SN, the streaming instability is 
more important. 



VI. SUMMARY 

In the paper above we characterized interaction 
of cosmic rays with balanced interstellar turbulence 
driven at a large scale. Our results can be summa- 
rized as follows: 

1. fast modes provide the dominant contribution 
to scattering of cosmic rays in different phases of in- 
terstellar medium provided that the turbulent energy 
is injected at large scales. This happens in spite of 
the fact that the fast modes are more subjected to 
damping compared to Alfven modes. 

2. As damping of fast modes depends on the angle 
between the magnetic field and the wave direction of 
propagation, we find that field wondering determined 
by Alfven modes affects the damping of fast modes. 
At small scales the anisotropy of fast mode damping 
makes gyroresonant scattering within slab approxi- 
mation applicable. At larger scales where damping 
is negligible the isotropic gyroresonant scattering ap- 
proximation is applicable. 

3. Transient time damping (TTD) provides an im- 
portant means of cosmic ray transport. Use of 5 func- 
tion resonance entails errors, and therefore, resonance 
broadening is essential for TTD. Our study shows that 
it is vital for low energy and large pitch angle scatter- 
ing. And it dominates scattering of cosmic rays in 
HIM and the partially ionized interstellar gas where 
fast modes are severely damped. 

4. Streaming instability is subjected to non-linear 
damping due to the interaction of the emerging mag- 
netic perturbations with the surrounding turbulence. 
The energy at which the streaming instability is sup- 
pressed depends on whether on the inverse cascade 
of magnetic energy as the instability gets more easily 
excited for low energy particles. 
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